(19) 



J) 



Europaisches Patentamt 
European Patent Office 
Office europeen des brevets 



(12) 



(43) Date of publication: 

02.05.2001 Bulletin 2001/18 

(21) Application number: 00309409.1 

(22) Date of filing: 25.10.2000 



III 

(ID EP 1 096 426 A1 

EUROPEAN PATENT APPLICATION 

(51) IntCI. 7 : G06T 11/00 



(84) Designated Contracting States: 


(72) 


Inventors: 


AT BE CH CY DE DK ES Fl FR GB GR IE IT LI LU 


• 


Besson, Guy W. 


MC NL PT SE 




Wauwatosa, Wisconsin 53213 (US) 


Designated Extension States: 


• 


Patch, Sarah Kathryn 


AL LT LV MK RO SI 




Saratoga Springs, New York 12866 (US) 


(30) Priority: 27.10.1999 US 428183 


(74) 


Representative: Pedder, James Cuthbert 






GE London Patent Operation, 


(71) Applicant: GE Medical Systems Global 




Essex House, 


Technology Company LLC 




12/13 Essex Street 


Waukesha, Wisconsin 53188 (US) 




London WC2R 3AA (GB) 



(54) Methods and apparatus for cone beam multi-slice CT correction 



(57) The present invention is, in one embodiment, a 
method for multi-slice, computed tomography imaging. 
The method includes steps of moving an (14) x-ray 
source through a trajectory; projecting an x-ray cone 
beam (16) from the moving x-ray source through an ob- 



ject towards a curved detector (20); and determining 
contributions of segments along the trajectory to voxels 
in a reconstruction volume of the object, including com- 
pensating the determined contributions for curvature of 
the detector and x-ray cone beam geometry. 
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Description 

[0001] This invention relates generally to methods and apparatus for reconstruction of imaging data, and more par- 
ticularly to cone beam correction for reconstruction of three-dimensional computed tomography imaging data. 

s [0002] In at least one known computed tomography (CT) imaging system configuration, an x-ray source projects a 
fan-shaped beam which is collimated to lie within an X-Y plane of a Cartesian coordinate system and generally referred 
to as the "imaging plane". The x-ray beam passes through the object being imaged, such as a patient. The beam, after 
being attenuated by the object, impinges upon an array of radiation detectors. The intensity of the attenuated beam 
radiation received at the detector array is dependent upon the attenuation of the x-ray beam by the object. Each detector 

10 element of the array produces a separate electrical signal that is a measurement of the beam attenuation at the detector 
location. The attenuation measurements from all the detectors are acquired separately to produce a transmission 
profile. 

[0003] In known third generation CT systems, the x-ray source and the detector array are rotated with a gantry within 
the imaging plane and around the object to be imaged so that the angle at which the x-ray beam intersects the object 

15 constantly changes. A group of x-ray attenuation measurements, i.e., projection data, from the detector array at one 
gantry angle is referred to as a "view". A "scan" of the object comprises a set of views made at different gantry angles, 
or view angles, during one revolution of the x-ray source and detector. In an axial scan, the projection data is processed 
to construct an image that corresponds to a two dimensional slice taken through the object. One method for recon- 
structing an image from a set of projection data is referred to in the art as the filtered back-projection technique. This 

20 process converts the attenuation measurements from a scan into integers called "CT numbers" or "Hounsfield units", 
which are used to control the brightness of a corresponding pixel on a cathode ray tube display. 
[0004] In at least one known multi-slice CT system, a "cone angle," or volumetric content of measured data, is very 
small. Therefore, this system currently processes 3-dimensional data using a 2-dimensional algorithm. By using the 
Feldkamp (FDK) algorithm, which is a simple perturbation of a 2-dimensional filtered back projection (FBP) algorithm 

25 for image reconstruction, excellent image quality is obtained for these relatively small cone angles. The FDK algorithm 
is not exact, however, and as the number of slices increases (for fixed slice thickness) cone angle and cone beam 
artifacts increase. 

[0005] It would be desirable to extend 2-dimensional CT fan-beam reconstruction algorithms to the cone-beam ge- 
ometry of a third-generation multi-slice CT imaging system. Such a reconstruction could be based upon a modified 

30 FDK algorithm, but the modifications would have to compensate for both a cylindrical (rather than planar) area detector, 
and a helical (rather than circular) source trajectory. Use of a curved detector array requires data interpolation along 
curved lines on the detector, the application of a new data filter, and of pre- and post-convolution weights. However, 
the helical source trajectory complicates voxel-driven back projection in that variable data (projection ray) redundancy 
conditions are encountered across a reconstructed image. For each voxel in a reconstruction volume, it is necessary 

35 to compute, for each source position, a ray passing from the x-ray source through the voxel. Thus, approaches must 
be found to address the handling of data redundancy and handling the variable number of rays that contribute to voxels. 
[0006] Two approaches to addressing the problems of data redundancy and the variable number of rays contributing 
to voxels are known. I n one approach , the helical pitch is limited so that at least two (interpolated) samples are obtained 
for each voxel. Extra data is thrown away, retaining only 2n worth of projections. The z-resolution available from con- 

40 jugate rays is ignored to keep the approach simple. However, in patient scanning, this, approach is impractical, because 
it severely limits patient coverage and results in increases radiation dose. A second approach improves on the first by 
providing better image quality (IQ) and reduced patient dose, while handling any practical pitch. However, a method 
implementing this approach to handling the above-mentioned problems requires that all data passing through a given 
voxel (for given source and fan angles) are simultaneously available. As a result, methods based on the second ap- 

45 p roach are impractical to implement. 

[0007] It would therefore be desirable to provide methods and apparatus that provide acceptable compromises be- 
tween improved image quality and practicality in third generation CT imaging systems. 

[0008] According to a first aspect of the invention, there is provided a method for multi-slice, computed tomography 
imaging, said method comprising: moving an x-ray source through a trajectory; projecting an x-ray cone beam from 
so the moving x-ray source through an object towards a curved detector; and determining contributions of segments along 
the trajectory to voxels in a reconstruction volume of the object, including compensating the determined contributions 
for curvature of the detector and x-ray cone beam geometry. 

[0009] Determining contributions of segments along the trajectory to voxels in a reconstruction volume of the object 
may comprise the step of filtering and back-projecting using a double integral over a fan angle and source angle for 
55 each point in the reconstruction volume. 

[0010] Back-projecting using a double integral may comprises back-projecting onto a polar reconstruction lattice. 
[001 1 ] The curved detector may be a cylindrical detector, and compensating the determined contributions for curva- 
ture of the detector and x-ray cone beam geometry may comprise the step of compensating the determined contribu- 
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tions for curvature of the cylindrical detector. 

[0012] Moving an x-ray source through a trajectory may comprise moving the x-ray source through a helical trajectory 
relative to the object. 

[0013] According to a second aspect of the invention, there is provided a method for multi-slice, computed tomo- 
5 graphic imaging, said method comprising: moving an x-ray source through a trajectory; projecting an x-ray cone beam 
from the moving x-ray source through an object towards a curved detector; and weighting, filtering, and back-projecting 
a signal representative of detected x-rays passing through the object to produce a tomographic image; the back-pro- 
jection being performed in a three-dimensional space, with compensation for x-ray cone-beam geometry. 
[0014] Back-projecting using a signal representative of detected x-rays may comprise the step of back-projecting 
10 onto a polar reconstruction lattice. 

[0015] The weighting may be performed prior to said filtering, and said filtering may be performed prior to said back- 
projection. 

[0016] The weighting, filtering, and back-projecting a signal representative of detected x-rays passing through the 
object to produce a tomographic image may comprise determining a reconstruction image function f(v) written as: 

15 

* L 2 (v,/7) (cos 2 y+sin 2 /cos 2 a) 



20 



25 



where: 

® represents convolution; 
S is a source to isocenter distance; 
30 v is a voxel with cylindrical coordinates (r, <p, k); 

= (v,p) is a z-elevation for a projection of v onto a detector; 
L is a voxel to source distance (in three-dimensional space); 

35 

g() denotes a fan-beam reconstruction convolution kernel; 
p* is data that has been interpolated onto a tilted plane; 
40 HSWQ denotes the half-scan weights; 

7 is a fan angle; 
p is a source angle; 

45 

a is an angle between a horizontal plane and the tilted planes; and 
r is a maximum fan angle of the fan beam. 
so [0017] The method may operate with the equation:- 



55 



g(/-r) = hijr-Y)^A{y)K(y-Y)h{y-y)B{y); 



where h{) is a parallel kernel and A, B, and Kare written as: 
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4(Y)= l±coglatenS, 
1 + tan y 



v 1 + cos z a tan 2 y 



r)-— 



tan 2 x 



to 



75 



20 



r-r 



K i y\ l + tan 2 (/-x) ^ . 

1 + cos 2 a tan 2 (/ - y ) ^sinja tan[cosa tan(/ -x)]} J 



and 

t|(y) = arctan[(cosa) x tan(y)]; 



and where r is a fan angle for a ray passing through a voxel being reconstructed. 

[0018] Moving the x-ray source through a trajectory may comprise the step of moving the x-ray source in a helical 
trajectory relative to the object. 
25 [0019] Projecting an x-ray cone beam from the moving x-ray source through an object towards a curved detector 
may comprise the step of projecting the x-ray cone beam from the moving x-ray source through an object towards a 
cylindrical detector. 

[0020] Weighting, filtering, and back-projecting a signal representative of detected x-rays passing through the object 
to produce a tomographic image may comprise determining a reconstruction image function /(v)=/(r, <J>, k) written as: 



30 



35 



40 



o -r L + K cos\a) 



where: 

B denotes an overscan angle; 
SW(p,y) denotes a member of a group consisting of overscan and underscan weight functions; 
S is a source to isocenter distance; 
45 v is a voxel with cylindrical coordinates (r, <p, k); 

% = % (v,p) is a z-elevation for a projection of v onto a detector; 
L is a voxel to source distance (in three-dimensional space); 

50 

g() denotes a fan-beam reconstruction convolution kernel; 
p* is data that has been interpolated onto a tilted plane; 
55 y is a fan angle; 

r is a fan angle for a ray passing through a voxel being reconstructed; 
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p is a source angle; - 

a is an angle between a horizontal plane and the tilted planes; and 
r is a maximum fan angle of the fan beam. 
[0021] The method may operate in accordance with the equation 



g{r-r) = 



\2 



Y-Y 



,sin(^)- 77(f)); 



h{y - y ) « A(r)K(r - y)h(y - y)b(y) ; 



where h{) is a parallel kernel and A, B, and Kare written as: 



20 



^ (Y)= 1 + cos 2 «tan 2 y . 
1 +tan y 



25 



f v 1 + cos 2 a tan 2 7 



tan 2 y 



30 



4r-r) = - 



l + tan 2 (/ -y) 
cos 2 atan 2 (y - y) 



\ 2 



Y-Y 



sinja tan[cosa tan(^ - / )]}> 



35 



and 



40 



45 



50 



55 



T|(y) = arctan[(cosa) x tan(y)]. 

[0022] The SW(p,y) may be an underscan weight function, and said weighting may be performed prior to said filtering, 
and said filtering may be performed prior to said back-projection. 

[0023] The SW(p,y) may be an overscan weight function, and said filtering may be performed prior to said weighting, 
and said weighting may be performed prior to said back-projection. 

[0024] The weighting may be performed prior to said filtering, and said filtering may be performed prior to said back- 
projecting. 

[0025] According to a third aspect of the invention, there is provided a multi-slice : computed tomography CT imaging 
system, said imaging system comprising a moveable x-ray source and a curved detector, said system configured to: 
move said x-ray source through a trajectory; project an x-ray cone beam from said moving x-ray source through an 
object towards said curved detector; and determine contributions of segments along the trajectory to voxels in a re- 
construction volume of the object, including compensating the determined contributions for curvature of said detector 
and x-ray cone beam geometry. 

[0026] The system may be configured to determine contributions of segments along the trajectory to voxels in a 
reconstruction volume of the object and the system may be configured to filter and back-project using a double integral 
over a fan angle and source angle for each point in the reconstruction volume. 

[0027] The system may be configured to back-project using a double integral comprises said system being configured 
to back-project onto a polar reconstruction lattice. 

[0028] The detector may be a cylindrical detector, and said system may be configured to compensate the determined 
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contributions for curvature of said detector and x-ray cone beam geometry and said system may be configured to 
compensate the determined contributions for curvature of said cylindrical detector. 

[0029] The system may be configured to move said x-ray source through a trajectory and said system may be con- 
figured to move said x-ray source through a helical trajectory relative to the object. 

5 [0030] According to a fourth aspect of the invention, there is provided a multi-slice, computed tomography CT imaging 
system, said imaging system comprising a moveable x-ray source and a curved detector, said system configured to: 
move an x-ray source through a trajectory; project an x-ray cone beam from said moving x-ray source through an object 
towards a curved detector; and weight, filter, and back-project a signal representative of detected x-rays passing 
through the object to produce a tomographic image; said system being configured to back-project in a three-dimensional 

10 space, with compensation for x-ray cone-beam geometry. 

[0031] According to a fifth aspect of the invention, there is provided a multi-slice, computed tomography CT imaging 
system in accordance with Claim 23 wherein said system being configured to back-project using a signal representative 
of detected x-rays comprises said system being configured to back-project onto a polar reconstruction lattice. 
[0032] The system may be configured to weight said signal prior to filtering said signal, and to filter said signal prior 

15 to said back-projecting. 

[0033] The system may be configured to weight, filter, and back-project a signal representative of detected x-rays 
passing through the object to produce a tomographic image and said system may be configured to determine a recon- 
struction image function f(v) written as: 



20 



25 



30 



* L (v, fi) (cos 2 y +sin 2 /cos 2 a ) 



where: 

® represents convolution; 
S is a source to isocenter distance; 
v is a voxel with cylindrical coordinates (r, <p, k); 
35 % = (v,p) is a z-elevation for a projection of v onto a detector; 

L is a voxel to source distance (in three-dimensional space); 
g() denotes a fan-beam reconstruction convolution kernel; 

40 

p* is data that has been interpolated onto a tilted plane; 
HSW() denotes the half-scan weights; 
45 y is a fan angle; 

p is a source angle; 

a is an angle between a horizontal plane and the tilted planes; and 

50 

r is a maximum fan angle of the fan beam. 
[0034] 27. A method in accordance with Claim 26 wherein 

55 
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g(y-r) = 



Y-r 



h{y - f ) * A{y)K{y - y)h(y - y )b(? ) ; 



10 



where is a parallel kernel and >4, B, and Kare written as: 



A(y)= l±cosloyanS ; 
1 + tan y 



75 



x 1 + cos 2 a tan 2 y 



*(f) = 



1 + tan 2 y 



20 



25 



and 



W ~\ l + tan 2 (y-f) 

M V - v ) = ; —, 

v 7 1 + cos- atan 2 (y-y) 



V 



r-r 



^sin|a tan[cosa tan(y - y )]} > 



3D 



tj(y) = arctan[(cosa) x tan(y)); 



and where r is a fan angle for a ray passing through a voxel being reconstructed. 

[0035] The system may be configured to move said x-ray source through a trajectory and said system may be con- 
figured to move said x-ray source in a helical trajectory relative to the object. 
35 [0036] The curved detector may be a cylindrical detector. 

[0037] The system may be configured to weight, filter, and back-project a signal representative of detected x-rays 
passing through the object to produce a tomographic image and the system may be configured to determine a recon- 
struction image function f(v)=/(r, k) written as: 



40 



45 



50 



55 



* J r L 2 +/r 2 cos(a) 



where: 

B denotes an overscan angle; 

SW(P,y) denotes a member of a group consisting of overscan and underscan weight functions; 

S is a source to isocenter distance; 

v is a voxel with cylindrical coordinates (r, <p, k); 

£ = £ (v,p) is a z-elevation for a projection of v onto a detector; 
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L is a voxel to source distance (in three-dimensional space); 
g() denotes a fan-beam reconstruction convolution kernel; 
p* is data that has been interpolated onto a tilted plane; 
7 is a fan angle; 

y is a fan angle for a ray passing through a voxel being reconstructed; 
p is a source angle; 

a is an angle between a horizontal plane and the tilted planes; and 
r is a maximum fan angle of the fan beam. 
[0038] The method may use the formula 



20 



g{r-r) = 



r-r 



h{ y - y ) « a{y)k{ y - Mr - y)b(y) ; 



25 



where h{) is a parallel kernel and A, B } and Kare written as: 



30 



A , v 1 + cos a tan y 

A(y) = 2 ; 

1 + tan y 



35 



/v I + cos 2 a tan 2 r 

b{y) = 



1 + tan 2 y 



40 



45 



; and 



K{r ~ r) ~ l + cos 2 atan 2 (y-y) 



K sin{a tan[cos a tan(y - y )]| , 



T|(y) = arctan[(cosoc) x tan(y)]. 



[0039] The SW(P,y) may be an underscan weight function, and said system may be configured to perform said weight- 
50 ing prior to said filtering, and said filtering prior to said back-projection. 

[0040] The SW(p,y) may be an overscan weight function, and said system may be configured to perform said filtering 
prior to said weighting, and said weighting prior to said back-projection. 

[0041] The system may be configured to perform said weighting prior to said filtering, and said filtering prior to said 
back-projecting. 

55 [0042] There is therefore provided, in one embodiment, a method for multi-slice, computed tomography imaging. 
The method includes steps of moving an x-ray source through a trajectory; projecting an x-ray cone beam from the 
moving x-ray source through an object towards a curved detector; and determining contributions of segments along 
the trajectory to voxels in a reconstruction volume of the object, including compensating the determined contributions 
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for curvature of the detector and x-ray cone beam geometry. 

[0043] The above described embodiment and others described in detail herein provide improved image quality in 
third generation CT imaging systems, even for cylindrical detectors and helical source trajectories. In addition, these 
embodiments can be implemented within the practicality constraints imposed by CT imaging hardware and patient 
5 dosage limitations. 

[0044] The invention will now be described in greater detail, by way of example, with reference to the drawings, in 
which:- 

[0045] Figure 1 is a pictorial view of a CT imaging system. 

[0046] Figure 2 is a block schematic diagram of the system illustrated in Figure 1 . 
to [0047] Figure 3 is a geometric representation of an embodiment of a third generation, four slice CT imaging system. 
[0048] Figure 4 is a geometric diagram showing conversion of FDK infinitesimals into third generation geometry. 
[0049] Figure 5 is a geometric representation of fan beam geometry. 

[0050] Figure 6 is a geometric representation of a third generation, cone-beam geometry with a cylindrical detector. 
[0051] Figure 7 is a graphical representation of a polar reconstruction lattice. 

'5 (0052) *■ igurc 8 is a geometric representation of a third generation geometry and source trajectory. 

[0053] Kciernng to Figures 1 and 2, a computed tomograph (CT) imaging system 10 is shown as including a gantry 
12 rcpr escntrfttve of a "third generation" CT scanner. Gantry 12 has an x-ray source 14 that projects a beam of x-rays 
16 towa-d h detector array 1 8 on the opposite side of gantry 12. Detector array 1 8 is formed by detector elements 20 
wheh icgc'.hc sense the projected x-rays that pass through an object 22, for example a medical patient. Detector 

20 array 1 8 may be labr icated in a single slice or multi-slice configuration. Each detector element 20 produces an electrical 
signal that represents the intensity of an impinging x-ray beam and hence the attenuation of the beam as it passes 
through patient 22. During a scan to acquire x-ray projection data, gantry 12 and the components mounted thereon 
rotate about a center of rotation 24. 

[0054] Rotation of gantry 12 and the operation of x-ray source 14 are governed by a control mechanism 26 of CT 
25 system 1 0 Control mochanism 26 includes an x-ray controller 28 that provides power and timing signals to x-ray source 
14 and a gantry motor controller 30 that controls the rotational speed and position of gantry 12. A data acquisition 
system (DAS) 32 in control mechanism 26 samples analog data from detector elements 20 and converts the data to 
digital signals for subsequent processing. An image reconstructor 34 receives sampled and digitized x-ray data from 
DAS 32 and performs high speed image reconstruction. The reconstructed image is applied as an input to a computer 
30 36 which stores the image in a mass storage device 38. 

[0055] Computer 36 also receives commands and scanning parameters from an operator via console 40 that has a 
keyboard. An associated cathode ray tube display 42 allows the operator to observe the reconstructed image and other 
data from computer 36. The operator supplied commands and parameters are used by computer 36 to provide control 
signals and information to DAS 32, x-ray controller 28 and gantry motor controller 30. In addition, computer 36 operates 
35 a table motor controller 44 which controls a motorized table 46 to position patient 22 in gantry 12. Particularly, table 
46 moves portions of patient 22 through gantry opening 48. 

[0056] In one embodiment of the present invention, a 2-dimensional CT fan-beam reconstruction method is extended 
to a cone-beam geometry on a third generation multi-slice CT A curved, cylindrical detector is used rather than a planar 
area detector, as is a helical, rather than a circular, source trajectory. The curved detector array produces data that is 

40 interpolated along curved lines of the detector by subsequent processing, and a data filter is used to filter the interpolated 
data using pre- and post-convolution weights. The helical source trajectory produces variable data (projection ray) 
redundancy conditions across the reconstructed image of the voxel-driven back-projection. Therefore, for each view 
in the reconstruction volume, a ray passing from the x-ray source through each voxel is computed. The geometry of 
an embodiment of a third generation, four-slice system is shown schematically in Figure 3. Gantry 12 rotates about a 

45 z-axis 50. Detector 18 surface 52 lies on a cylinder, with four rows, each having 888 detector elements 20. Row data 
represent rays 54 that do not lie in horizontal reconstruction planes. A "cone angle" 56 with respect to a gantry plane 
(which itself is defined as a plane orthogonal to a z-axis that contains the source 14 point and bisects the detector 18 
cylinder) represents a degree of inconsistency in cone beam data as compared to two-dimensional CT data. Planes 
containing a focal spot or point 58 intersect detector surface 52 along arcs 60 that drop off as cos(y). A ray 62 that 

50 passes through voxel V1 hits detector array 18, so that data representing voxel V1 is estimated by interpolation. A ray 
64 through voxel V2. however, does not hit detector array 18, so data representing voxel V2 is estimated by extrapo- 
lation. 

[0057] In one embodiment, a compromise is effected between improved image quality and practicality by combining 
voxel-driven back-projection, data extrapolation, and weighting methods, such as half-scan, underscan, and overscan, 
55 to handle data redundancy conditions and smooth out discontinuities at the 0, 2n source angle interface for helical 
trajectories. Use of extrapolation allows high scanning pitches to be used by providing a constant number of projection 
measurements for all voxels for a given source position and by ensuring that no measurement is discarded. Extrapo- 
lation thus eliminates cone-angle and detector size-related variable voxel sampling of source positions. Source-related 
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data redundancies and associated vanishing weights (half-scan weights in high spec -MS) pitches, overscan or un- 
derscan in high quality (HQ) pitches) are used to manage data inconsistencies result! om the combination of cone- 
angle and helical scanning. To enable independent projection processing, HS and HC „ (ches are calculated such that 
at least one ray or two rays, respectively, contribute to a given voxel. 
5 [0058] In one embodiment, a Feldkamp (FDK) reconstruction algorithm is modified for use in third-generation multi- 
slice CT imaging systems. FDK is a perturbation of 2-dimensional filtered back-projection (FBP), which calculates a 
double integral over a fan angle, 7, and a source angle, p, for each point in a reconstructed image. In FDK, a contribution 
of each infinitesimal segment along a source trajectory to a voxel in a reconstruction volume is computed as if the 
segment were contributing to the 2-dimensional reconstruction in a plane defined by the infinitesimal segment and the 
10 point/voxel. Therefore, FDK is derived from 2-dimensional FBP by a change of variables. Referring to Figure 4 and 
voxel V, S is a source 14 to isocenter 82 distance; S + D is a source 14 to detector 18 distance; 0 is a source angle, 
P' is a source angle in a tilted plane; y is a fan angle; yMs a fan angle in the tilted plane; and a is an angle between a 
horizontal and a tilted plane, r is a radial vector, and r* is a radial vector in the tilted plane under consideration (i.e., 
tilted with respect to the gantry plane.) Embodiments in which flat panel x-ray detectors are used do not require inter- 
ns polation onto cos(y) arcs; i.e., filtering along detector rows is equivalent to filtering co-planar data for such panels. 
Therefore, the changes of variables differ somewhat in embodiments using a flat panel x-ray detector 
[0059] After interpolating projection data onto planar sets, projection rays are distributed within the fan according to 
a fan-beam parameterization written as: 

20 

s* dp* =s dp 

tan/' = tan/ cosa 

25 

where y 1 represents the fan angle in the tilted plane, and dp 1 represents the infinitesimal arc length in the tilted plane 
indicated by r 1 as shown in Figure 4. Thus, 

t . t cosa dy /0 v 
30 cosy dy = — * ^ (2) 

(cos y + cos a cos y) 

[0060] Plugging (1) and (2) into the 2-dimensional FBP reconstruction formula thus results in an embodiment that 
exactly implements FDK on a third generation CT imaging system. 
35 [0061] An approximation leads to an embodiment of a shift-invariant filter that very closely approximates the true, 
shift-variant filter, but which is simpler and more readily accommodates convolution in the Fourier domain. One em- 
bodiment of a reconstruction algorithm of the present invention utilizing an exact shift-variant filter for FDK and its shift- 
invariant approximation is described as follows. 

[0062] A third-generation FDK fan-beam parameterization in tilted planes is given by (3): 



40 



45 



50 



r^y') = arctan[(cos#) x tan(/')] . (3) 
[0063] A shift-invariant approximation to a shift-variant kernel (written omitting the superscript t) is written: 



55 [0064] Accordingly, the following equation is written by substituting fort| in equation (4) above: 
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s{r-r) = 



v 



r-y 



h[ r - f ) « A( r )K{r - y)h{ r - rMr) (5) 



where /?() is a parallel kernel and A, B, and K are written as: 



10 



1 +cos 2 a tan 2 y . 



pre-convolution weight: A(y) = 2 

1 + tan y 



(5a) 



75 



post-convolution weight: £(y ) = 



_ x 1 + cos 2 a tan 2 ^ 



1 + tan 2 y 



(5b) 



20 [0065] The kernel terms that multiply the expression for the parallel kernel h are written as: 



r-y 



25 



kL y)= L!L^!kzZ) _. 

1 + cos z a tan* (y - y) ^ s in{a tan[cosa tan(y - y )]}, 



(6) 



[0066] By replacing the shift-variant kernel g by the shift-invariant approximation above into equation (8), an exact 
30 third-generation FDK algorithm embodiment is obtained. 

[0067] In another embodiment useful for small cone-angles, such as those associated with current multi-slice scan- 
ners, an approximation that directly renders the filtering shift-invariant (and therefore enables convolution) is utilized. 
This embodiment is useful for high speed mode. For small angles, the convolution is simplified by approximating fan- 
beam parameterization by: 

35 



y* = y cosct 



(7) 



40 



where oc= a(v, P). The reconstruction algorithm is then written as: 



45 



o L ( v > P) (cos 2 ^+sm 2 /cos* a ) 



or, in a half-scan weighted embodiment, 



50 



55 



^ L (v y P) (cos 2 y+sur ycos 2 a J 



where: 
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® represents convolution; 



r is a maximum fan angle of the fan beam; 



5 



S is a source to isocenter distance; 



v is a voxel with cylindrical coordinates (r, <p, k); 



% = s (v p) is a z-elevation for a projection of v onto a detector; 



10 



L is a voxel to source distance (in three-dimensional space); 



px) denotes a fan-beam reconstruction convolution kernel; 



15 



r is data that has been interpolated onto a tilted plane; and 



HSW denotes the half-scan weights. 

[0068] Tc process image data in accordance with equation (9), x-ray source 14 is moved through a trajectory, and 
20 a cone beam 1 6 is projected from x-ray source 1 4 through an object or patient 22 towards curved detector 18. A signal 
representative of detected x-rays 16 passing through object 22 is processed by Imaging system 10, for example, by 
DAS 32 imaging reconstructor 34, and computer 36, to produce a tomographic image that appears on CRT display 
42. (As used herein, the "signal representative of detected x-rays" is intended to be broadly construed. For example, 
embodiments in which parallel data lines for detector elements 20 are used to transmit the signal from detector 18 to 
25 DAS 32 arc mtondod to be included, as are embodiments which multiplex data on a single line.) The processing of the 
signal in accordance with equation 9 comprises applying half scan weights to pre-processed data, priorto reconstruction 
filtering These weights account for the fact that several line integrals have been measured twice, while most have 
been measured only once. Next, reconstruction filtering is performed, as represented by the convolution of g and 
p ■ HSW tn equation 9. Next, back-projection is applied in three-dimensional space, with compensation for x-ray cone 
30 beam geometry, by summing contributions from each projection to a given pixel. This back-projection is represented 
by the integral sum in equation 9. 

[0069] In an embodiment useful for high quality modes, interpolation of projection data along cos(y) arcs is eliminated. 
Accordingly, projection data are filtered along detector rows. The parameterization approximation (8) is retained. 
[0070] A derivation of one embodiment of an FDK algorithm of the present invention approximating FDK for third 

35 generation CT imaging systems 1 0 without cos(y) detector interpolation involves two approximations. The approxima- 
tions are first that the cone angle is small, so that convolution filtering is possible, and second, that the data are filtered 
along detector rows (i.e., without interpolation along cos(y) arcs). Figure 5 is a geometric representation of fan-beam 
geometry. All angles and distances refer to a central plane of gantry 12, orthogonal to a z-axis. A voxel M is shown 
projected onto the gantry 12 central plane. Considering reconstruction of voxel M with a multi-slice detector 18, a 

40 change of geometry from fan-beam to cone-beam results in a source to isocenter distance S changing into a source 
to z-axis S' given by: 



where £ is measured on the z-axis, and 

the source to detector distance E= S + D (as measured on a ray through the z-axis) is changed into £* given by: 



45 




50 



E' 2 =[£ 2 +£' 2 ]. 



where £ is measured on detector 1 8. 
55 [0071] Starting from a 2D reconstruction equation for fan-beam data, and in the notation of Figure 5: 
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/M= J4-J p(frr)g(r-r)cos r dy dp (10) 



-r 



with: 

10 



where h(u) denotes a parallel -geometry kernel. 
15 [0072] Using a multi-slice detector, the geometry changes from fan-beam to cone-beam in reconstruction of a voxel 
M. As a result, a source to isocenter distance S is changed into the source to z-axis distance S* written as: 



20 



where £ is measured on the z-axis. 

[0073] Also, the source to detector distance (as measured on a ray through the z-axis), is changed from E to Et 
written as: 

25 



E' 2 =[£ J +^ 2 ]. 



30 where is measured on the detector. 

[0074] Figure 6 is a geometric representation of a third-generation, cone-beam geometry with a cylindrical detector 
18. A ray 66 with a cone-angle a with respect to a plane of gantry 12 is shown, the plane being defined by O, x, andy 
in Figure 6. Also shown is an associated ray 68 (for the same elevation on detector 18) passing through an axis of 
rotation z. A rotation angle dp' such that a small rotation dp around z is equivalent to dp' around r is written as: 

35 

dp # --^=£fl>. „ (11) 

40 [0075] A fan-angle y in a plane of the gantry (orthogonal to z) is now replaced by a fan-angle y* in a tilted surface: 

y*= ~ y = 7 cos (a), 
E 

45 

(exactly under the assumptions above). 
[0076] Accordingly: 

55 [0077] Voxel-dependent geometric weight factor is determined. As data are filtered along detector rows, for a given 
voxel M, an associated £ to be used in the calculation of dp' above is found. Given 
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5 it is thus written: 

^ = |k. (13) 
10 [0078] Putting all the terms together, an equation for a third-generation FDK reconstruction algorithm is written: 



/(r.#.*)-7|7^T^-yy^.rtco{|- r) dy dfi (14) 



[0079] Alternatively, an equation for an embodiment of a modified FDK algorithm is written as: 

20 

/M'Tl ri 5 - \ , g{r-r)pV*r,?h4r cosfe)) «f y (15) 

oil +/c" cos(a) V J 



or, in an overscan or underscan weighted embodiment, 

30 

/(r,*.*)=T J 2 5 /^^ (r-yK Og^^OcosCrcosCa)^ (16) 

35 

where B denotes the overscan angle, and SW(0,-y) denotes a weight function that is either an overscan weight function 
or an underscan weight function. (For underscan, B=0). 

[0080] To process image data in accordance with equation (16), x-ray source 14 is moved through a trajectory. A 
cone beam 1 6 is projected from x-ray source 1 4 through an object or patient 22 towards curved detector 1 8. A signal 

40 representative of detected x-rays 16 passing through object 22 is processed by imaging system 10 by, for example, 
DAS 32, imaging reconstructor 34, and computer 32, to produce a tomographic image that is displayed on CRT display 
42. The processing of the signal in accordance with equation (16), when SW(fi,i) is an underscan weighting function, 
comprises weighting, filtering, and back-projecting the data, in that order. Back-projection is performed in three-dimen- 
sional space, with compensation for x-ray cone beam geometry. When StV(p,y) is an overscan weighting function, the 

45 weights depend only upon a source angle. In this case, the weights account for the fact that projection data have been 
acquired over a range of source angles larger than 27t, and accordingly, several line Integrals have been measured 
three times. Thus, filtering is performed before weighting, and weighting is performed before back-projecting. 
[0081] In one embodiment, an FDK algorithm is adapted to helical trajectories. In known multi-slice CT systems, 
most imaging is helical. Helical scanning is possible by using a modification in which all z distances are calculated with 

so respect to a gantry central plane. However, appropriate attention must be given to obtaining an appropriate number of 
data samples for each voxel in the reconstructed image and handling of data redundancy. Providing at least two samples 
from the detector using only interpolation leads to very slow pitch and 2 to 3 samples per voxel; on the other hand, 
handling sampling discrepancy by simply throwing away a third sample is not acceptable from a patient dose standpoint. 
Therefore, in one embodiment, extrapolation is used to enable faster pitches while taking into account all projection 

55 samples. Standard weighting techniques are used to handle data redundancy. Specifically, a HQ mode utilizing this 
embodiment accounts for two or more samples per voxel. Data redundancy from a.third sample is handled by employing 
overscan or underscan weights to smooth away data discrepancies resulting from the helical trajectory. In an HS mode 
embodiment, one or two samples per voxel are used. Data redundancies associated with the second sample are 
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handled using half-scan weights to feather data inconsistencies resulting from the helical trajectory and cone-angle. 
[0082] In one embodiment, for an eight-slice system, a 7:1 pitch is used to provide a best cone-artifact abatement 
measure. A modified FDK algorithm given by equation (15) is used with an additional modification making the gantry 
plane source-angle dependent. Also, the back-projection was extended from 2rcto (S/7+e)x2n using overscan weights, 
5 where e is a parameter of the method. The use of overscan weights avoids data inconsistencies of helical algorithm 
implementations at source angles 0 and 2n relating to a z-motion of the patient with respect to the source. (Imaging 
without handling this inconsistency results in streak artifacts in a direction of the source position at angle 0.) Useful 
overscan functions include those written as: 



where x varies between 0 and 1 in the interval considered, and 5 is a parameter. 

[0083] In one embodiment, high pitch reconstructions are implemented by applying !4 scan weights to the modified 
25 FDK algorithm described above, equation (8). In another embodiment, the data are back-projected on a three-dimen- 
sional Cartesian grid. 

[0084] In one embodiment of the present invention, x-ray source 14 is moved through a helical trajectory relative to 
object 22. An x-ray cone beam 16 is projected from source 14 through object 22 to a curved, cylindrical detector 18. 
Processing of detector data to produce CT images in imaging system 1 0 is provided by DAS 32, imaging reconstructor 
30 34, and computer 36, for example, and an image is displayed on CRT display 42. This processing includes determining 
contributions of segments along the trajectory to voxels in a reconstruction volume of the object, including compensating 
the determined contributions for curvature of the detector and x-ray cone beam geometry. This determination comprises 
filtering and back-projecting using a double integral over a fan angle and source angle for each point in the reconstruc- 
tion volume. 

35 [0085] Referring to the graphical representation of a polar reconstruction lattice shown in Figure 7, an exemplary 
embodiment of the present invention was implemented using a reconstruction done on a polar lattice using MATLAB. 
The dimensions of the axes of Figure 7 are in cm. for the exemplary embodiment. Forty-one of 984 source positions 
70 are shown. An undersampled reconstruction lattice 72 is also shown. The image was interpolated onto a Cartesian 
lattice using MATLAB's cubic interpolation call. A polar lattice is an appropriate geometry for a third generation system 

40 having a circular source trajectory. An angular lattice increment equal to the angular increment of CT scanner 10 
between views was selected, d9 = dp = 2iE/N views . so tnat dr~dx, where dx is the size of a pixel of the final image. 
Thus, a distance between reconstruction lattice points is roughly dx at outer edges of the FOV, decreasing to zero at 
the center. The polar lattice simplified programming, and further permitted precomputation of all trigonometric functions 
used in the back-projection. 

45 [0086] In Figure 7, R is a source trajectory radius and kdQ represents a source angle. A pixel under consideration is 
at a radial distance r and an angle jdQ. With this in mind, it is relatively simple to compute L, a source-pixel distance. 
The only trigonometric functions needed to calculate these distances L are {cos(i otf)} for M>, ... (number of views - 1). 
An extension to a 3d back-projection is also simple, and is given by L 3d (Z- 2 +cte 2 ) 1/2 . Pre-computation significantly re- 
duced the algorithm's floating point operation (FLOP) count, and efficiency could also be improved by avoiding back- 

50 projection onto every lattice point near the origin. (Selective back-projection onto some points in the polar lattice in- 
creases the programming burden.) Referring to the geometric representation of third generation and source trajectory 
shown in Figure 8, 0 th order extrapolation was used when the required data in the voxel-driven back-projection fell 
above or below the 8-slice detector. First order or higher order, or non-polynomial methods, might be used for the 
extrapolation. A source trajectory 74 and detector surface 52 are shown. A region of reconstruction is shown at 76. 

55 Half-scan data was collected only between view angles -(tc/2 + T) and +(n/2 + T). A ray 78 through a point 80 at view 
angle -(n/2 + T) does not hit detector 18 surface 52 for helical pitches above N^ lce J2, and hence, data for ray 78 is 
extrapolated. However, a conjugate ray 82, collected a view angle 0, is always measured by detector 18 with cone 
angle = 0. In one embodiment illustrated in Figure 8, R= 54.1 cm; D = 40.8 cm; FOV =48 cm; r = 27.4°; = pitch 



10 



f(x) = 3x 2 -2x 3 , 



where x varies between 0 and 1 in the interval considered; 



15 




20 
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x slice width x (tc/2 + y)/(2n); and k = ray height on detector = -z^^ x {R + D)/L^ 1 .06 x pitch x slice width. 
[0087] From the preceding description of various embodiments of the present invention, it is evident that the embod- 
iments described herein lead to simple changes in data processing that makes reconstruction more robust to large 
cone angles. Data are filtered "on the fly," with only minor changes to standard reconstruction chains. The embodiments 
5 described herein should be contrasted with previously known adaptations of the FDK algorithm to third-generation CT 
with cylindrical detectors that fail to account for geometry and thus filter non-coplanar data. In addition, the embodiments 
described herein do not require simultaneous data availability for ail rays passing through a voxel for a given 2D Radon 
point as does at least one other known reconstruction method. The embodiments described herein also fully use patient 
dose while providing increased patient coverage. Artifacts are effectively eliminated by the use of over-scan weights 
10 or underscan weights (respectively half-scan weights), which, in combination with the use of extrapolation, makes the 
use of high pitches practical in all imaging modes without wasting patient dose. The embodiments described herein 
use relatively few views as compared to at least one previously known four slice helical reconstruction algorithm, and 
accordingly lead to fast reconstruction and improved temporal resolution. It will be recognized by those skilled in the 
art that the method and apparatus embodiments of the present invention described herein provide acceptable corn- 
's promises between Improved image quality and practicality in third generation CT imaging systems. Moreover, these 
embodiments can be implemented within the practicality constraints imposed by CT imaging hardware and patient 
dosage limitations. 

[0088] Although particular embodiments of the invention have been described and illustrated in detail, it is to be 
clearly understood that the same is intended by way of illustration and example only and is not to be taken by way of 

20 limitation. For example, embodiments of the invention are applicable to 8. 1 6, or other numbers of detector rows. Also, 
the embodiments described herein use detector row-to-row interpolation for a given source position. However, further 
improvement in image quality at small cone angles can be obtained, with somewhat increased implementation com- 
plexity, by using conjugate ray interpolation/extrapolation to zero-cone angle. In addition, the CT system described 
herein is a "third generation" system in which both the x-ray source and detector rotate with the gantry. Many other CT 

25 systems including 'fourth generation" systems wherein the detector is a full-ring stationary detector and only the x-ray 
source rotates with the gantry, may be used if individual detector elements are corrected to provide substantially uniform 
responses to a given x-ray beam. 

30 Claims 

1 . A method for multi-slice, computed tomography imaging, said method comprising: 

moving an x-ray source (14) through a trajectory (74); 

35 

projecting an x-ray cone beam (1 6) from the moving x-ray source through an object towards a curved detector 
(20); and 

determining contributions of segments along the trajectory to voxels in a reconstruction volume (76) of the 
40 object, including compensating the determined contributions for curvature of the detector and x-ray cone beam 

geometry. 

2. A method in accordance with Claim 1 and further wherein determining contributions of segments along the trajec- 
tory (74) to voxels in a reconstruction volume (76) of the object comprises the step of filtering and back-projecting 

45 using a double integral over a fan angle and source angle for each point in the reconstruction volume. 

3. A method in accordance with Claim 2 wherein back-projecting using a double integral comprises back-projecting 
onto a polar reconstruction lattice. 

50 4. a method for multi-slice, computed tomographic imaging, said method comprising: 

moving an x-ray source (14) through a trajectory (74); 

projecting an x-ray cone beam (1 6) from the moving x-ray source through an object towards a curved detector 
55 (20); and 

weighting, filtering, and back-projecting a signal representative of detected x-rays passing through the object 
to produce a tomographic image; 
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the back-projection being performed in a three-dimensional space, with compensation for x-ray cone-beam 
geometry. 

5. A method in accordance with Claim 4 wherein back-projecting using a signal representative of detected x-rays 
5 comprises the step of back-projecting onto a polar reconstruction lattice (72). 

6. A method in accordance with Claim 4 or 5, wherein said weighting is performed prior to said filtering, and said 
filtering is performed prior to said back-projection. 

10 7. A multi-slice, computed tomography CT imaging system (1 0), said imaging system comprising a moveable x-ray 
source (1 4) and a curved detector (20), said system configured to: 

move said x-ray source through a trajectory (74); 

15 project an x-ray cone beam (1 6) from said moving x-ray source through an object towards said curved detector; 

and 

determine contributions of segments along the trajectory to voxels in a reconstruction volume (76) of the object, 
including compensating the determined contributions for curvature of said detector and x-ray cone beam ge- 
20 ometry. 

8. A multi-slice, computed tomography CT imaging system (1 0) in accordance with Claim 7 and further wherein said 
system being configured to determine contributions of segments (20) along the trajectory (74) to voxels in a re- 
construction volume (76) of the object comprises said system being configured to filter and backproject using a 

25 double integral over a fan angle and source angle for each point in the reconstruction volume. 

9. A multi-slice, computed tomography CT imaging system (10), said imaging system comprising a moveable x-ray 
source (14) and a curved detector (20), said system being configured to: 

30 move an x-ray source through a trajectory (74); 

project an x-ray cone beam (16) from said moving x-ray source through an object towards a curved detector; 
and 

35 weight, filter, and back-project a signal representative of detected x-rays passing through the object to produce 

a tomographic image; 

said system being configured to back-project in a three-dimensional space, with compensation for x-ray cone- 
beam geometry. 



40 



45 



50 



55 



10. A multi-slice, computed tomography CT imaging system (10) in accordance with Claim 9 wherein said system 
being configured to back-project using a signal representative of detected x-rays comprises said system being 
configured to back-project onto a polar reconstruction lattice (72). 
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